Optimal control of entangling operations for trapped ion quantum computing 
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Optimal control techniques are applied for the decomposition of unitary quantum operations into 
a sequence of single-qubit gates and entangling operations. To this end, we modify a gradient-ascent 
algorithm developed for systems of coupled nuclear spins in molecules to make it suitable for trapped 
ion quantum computing. We decompose unitary operations into entangling gates that are based 
on a nonhnear collective spin operator and complemented by global spin fiip and local light shift 
gates. Among others, we provide explicit decompositions of controUed-NOT and Toffoli gates, and 
a simple quantum error correction protocol. 



PACS numbers: 03.67-a, 32.80.Qk, 37.10.Ty 



I. INTRODUCTION 

Choosing the laws of quantum physics as the physi- 
cal basis for constructing models of computation 1] al- 
lows for solving certain computational problems more ef- 
ficiently as in models based on classical physics 0. In 
the quantum circuit model, information is encoded in 
quantum bits (qubits) and manipulated by applying uni- 
tary operations acting on the joint state space of the 
qubits. It has been shown that arbitrary unitary opera- 
tions can be broken down into sequences of elementary 
gate operations, consisting of single-qubit operations and 
entangling operations acting on pairs of qubits Q . It is a 
non-trivial task to find optimum decompositions of uni- 
tary operations into a minimum number of elementary 
gates. Often, the controUed-NOT (CNOT) gate opera- 
tion is chosen as the entangling operation. However, it 
has been shown that almost any entangling operation can 
be used for this purpose as well 4]. 

In experiments processing quantum information, the 
available physical processes determine the choice of the 
entangling gate. In the case of trapped ions manipulated 
by coherent laser light 0], a qubit is realized by encod- 
ing quantum information in a pair of long-lived internal 
states, consisting of hypcrfine or Zeeman ground states, 
or in a combination of a ground state and a metastable 
state with an energy difference of a few electron volts. 
Gates acting on a single qubit are achieved by lasers cou- 
pling the qubit states by either dipole-forbidden single- 
photon transitions or Raman transitions. Pairs of ions 
are entangled by qubit-qubit interactions mediated by a 
coupling to a vibrational mode of the ions' motion in the 
trap. 

To achieve a gate on either a specific ion or a pair of 
ions within an ion string, two strategies are pursued: 

1. The laser beam is tightly focussed so that it inter- 
acts only with a single ion at a time @. Then, a 
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CNOT gate or an equivalent entangling gate be- 
tween an arbitrary pair of ions is achieved by a se- 
quence of pulses with the laser addressed to either 
one or the other ion of the pair Q . 

2. Alternatively, a laser with wider beam diameter is 
employed in combination with ions held in a seg- 
mented ion trap. In this approach, trap potentials 
are dynamically transformed to enable the trans- 
port of one or two ions into the interaction region 
with the laser beam. To entangle a pair of ions i, j, 
a bichromatic laser field is used to realize either a 
conditional phase gate [8] induced by an effective 

Hamiltonian Hpo cx cri^Vi"'-' or a M0lmer-S0rensen 
gate induced by Hms « Vi"'"' where ait^ de- 
notes a Pauli spin operator cr ■ n acting on the /c'th 



Recently, the latter interaction has been used to entangle 
a pair of ions with fidelities of up to 99.3(1)%, with a cou- 
pling mediated by the longitudinal center-of-mass mode 
of the ion string 10]. Employing the same kind of inter- 
action to iV > 2 ions would realize a unitary operation 

^Msi^) = expi-i^S^"^) with = J2k=i '^x^ yielding an 
equal pairwise coupling between all ions in the string. 
In addition, using a single frequency resonant with the 
qubit transition instead of a bichromatic field coupling to 
the transitions' motional sidebands, the same laser beam 
could induce spin flips on all ions described by the uni- 
tary Ux{9) — exp(— il^a;). This opens up the interest- 
ing prospect of performing arbitrary unitary operations 
by complementing these unitaries by single-qubit phase 
shift gates ui^\9) = exp(— i|tTi'^^) that could be induced 
by a tightly focussed off-resonant laser beam interacting 
only with the /c'th ion. In this approach, no interferomet- 
ric stability is required between the optical path lengths 
of the focussed and the wide beam. Moreover, the use of 
light shift gates facilitates addressing a single qubit in a 
string of ion qubits without introducing unwanted state 
transformations on the neighboring ions since the phase 
shift is proportional to the intensity of the laser field as 
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compared to single-qubit spin flip gates where the rota- 
tion angle 9 is proportional to the field amplitude . In 
this way, unitary operations could be realized on a small 
group of ions without the need to split and rearrange the 
ion string in-between ^12,] and with more modest require- 
ments regarding the spatial mode profile of the tightly 
focussed laser beam. 

In this paper, we will use optimal control techniques 
to find decompositions of 7V-qubit gates into unitary op- 
erations induced by the set of Hamiltonians 

5 = {5.^5.,aW,ai2),...,af)}. (1) 

In contrast to similar applications of optimal control in 
nuclear magnetic resonance (NMR) experiments [isj and 
in laser-induced femto-chemistry |14, 15], we are inter- 
ested in the case where only one of the Hamiltonians is 
applied at a time. As a consequence, the optimal control 
algorithm is required to find decompositions of a given 
quantum gate by sequences of laser pulses that interact 
either with all ions in the same way or with an individual 
ion. 

The goal of this paper is to find gate decompositions 
of interest for current state-of-the-art ion trap experi- 
ments that arc more efficient than decompositions based 
on gates acting on only one or two qubits at a time. 



II. BASIS SET OF OPERATIONS 

The mean-field interaction Sx'^ acting on a string of N 
ions entangles each ion qubit with each other ion qubit. 
Thus, on the one hand, application of this Hamiltonian 
endows us with the power to entangle arbitrary pairs of 
qubits. On the other hand, it induces entangling inter- 
actions that are not always desired. The situation we 
are encountering is somewhat similar to the one in NMR 
quantum computing where the system Hamiltonian Hgys 
consisting of spin-spin interactions and chemical shifts 
leads to a system dynamics that needs to be controlled 
by radio-frequency fields interacting with single spins at 
a time. For this purpose, techniques have been devel- 
oped for selectively switching off certain spin-spin inter- 
actions by refocussing techniques [iGj . In the scenario 
we envision for the ion trap system, the role of the radio- 
frequency fields is taken over by laser pulses inducing 
single-qubit gate operations that are intermittently ap- 
plied to particular ions. 

Refocussing techniques are also applicable to the 
trapped ion system. For example, to entangle qubits 
1 and 2 in a system of three qubits, the entangling 
pulse could be split into two parts and interleaved 
with a refocussing pulse to obtain the sequence U = 

U^siTT/A)ui^\7r)U§s{T:/4). The light shift pulse on 
the third qubit flips its phase and effectively reverses the 
entangling interactions with qubits 1 and 2. Substitut- 
ing the last pulse by its inverse, results in a sequence 
U = [/j^5(-7r/4)l7p^(7r)t/j^5(7r/4) that entangles qubits 



1 and 2 with qubit 3 without inducing entangling inter- 
actions between 1 and 2. 

The set of Hamiltonians (H]) is sufficient for gen- 
erating arbitrary unitary operations on a string of 
N qubits as can be shown by an explicit construc- 
tion: A single-qubit x-rotation acting on qubit k 
is generated by the spin echo sequence ui^\9) — 
ui'"\-TT)Ux{-e/2)ui''\TT)Uxie/2). Arbitrary single- 
qubit gates can then be performed when combining this 
operation with single-qubit phase shift gates ui^\9). 
Similarly, an operation f/j^g^*" that entangles all ions 
with each other except for qubit k is produced from the 
A^-qubit entanghng gate U^g{9) by the pulse sequence 
Um7' = ui''\-7r)U^,si(^/2)U^.''\n)U^9/2). Sub- 
stituting U^jg in the above sequence by C^ms^*"' ^ simi- 
lar sequence is obtained that entangles all qubits except 
for two, and by induction, a two-qubit entangling gate is 
constructed between any pair of qubits m and n which to- 
gether with arbitrary single qubit gates forms a universal 
set of gates. 

While this construction shows that in principle arbi- 
trary unitary operations are realizable by pulse sequences 
generated from the set S, it is of no practical use. For 
the implementation of A^-qubit gate operations, we are 
interested in finding pulse sequences that minimize gate 
errors occurring in ion trap quantum computing. There- 
fore, we will be searching for sequences having either a 
minimum number of (entangling) pulses or a minimum 
length in terms of the sum of pulse angles 9n of the indi- 
vidual pulses. 



III. OPTIMAL CONTROL OF UNITARY 
TRANSFORMATIONS 

Optimal control techniques have been applied to the 
problem of generating specific unitary transformations 
[Tsl [Tsl . [l7j with applications to systems as differ- 
ent as NMR, neutral atoms in optical lattices, Josephson 
junction qubits and trapped ions [H, [H, [l^, [23| • To find 
decompositions of entangling ion trap gates in terms of 
pulses generated by Hamiltonians Hk from S, we mod- 
ify a gradient-ascent algorithm that was developed by 
Khaneja et al. [11] in the context of NMR experiments. 
In their approach, a unitary transformation C/targct was 
searched for by constructing a unitary operation 

M M 

U=l[U„,= Y[ (2) 

m— 1 m— 1 

where Hk € 5, that maximized the performance function 
^{{ukm}) = lTr(;7'^C/targct)l^- The authors uotcd that for 
small time increments At the calculation of the gradient 

m 

(9$ iAt 

« -2Re{Ti{ — W^HkV^)TT{W^Vmy), (3) 
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with Wm = UI^^^^^Un ■ ■ - Um+i and Vm = Um---Ui, 
could be efficiently carried out requiring only about 3M 
matrix multiplications and about KM calculation of 
traces. Then, a gradient-based algorithm was devised to 
increase the value of the performance function by modi- 
fying the control amplitudes 

Ukm Ukm + (4) 

OUkm 

using a suitable step size e. Repeated application of 
the gradient calculation followed by updating the con- 
trol amplitudes maximized the performance function and 
resulted in a unitary transformation realizing the target 
operation ?7targot- 

There are a few important differences between coupled 
spin systems in NMR and in laser-manipulated strings 
of trapped ions. In the NMR context of ref. [3], the 
product of unitaries in ^ arises from a discretization of 
the time variable that required the approximation of a 
continuous control amplitude Uk {t) by a stepwise contin- 
uous function with values Ukj. Restrictions to the values 
of Uk{t) are only due to technical requirements like am- 
plitude or bandwidth limitations of the radio-frequency 
equipment used for producing the control fields. In ion 
trap experiments, however, limitations exist for the val- 
ues that the functions Uk{t) can take on because the si- 
multaneous application of different control Hamiltonians 
is either technically very challenging or physically impos- 
sible. In the former case, the simultaneous application of 
single qubit phase shift gates to more than a single ion 
would require control of the spatial profile of a laser beam 
inducing phase shifts by the ac-Stark effect. More impor- 
tantly, in the latter case, the entangling interaction Sx"^ is 
produced by an effective Hamiltonian that precludes the 
simultaneous application of single-qubit phase shifts. As 
a consequence, the control functions need to satisfy the 
condition Uk{t)ui{t) ~ for all fc, I at all times t. There- 
fore, the unitary transformation is naturally decomposed 
into a product of unitaries and Eq. ([2|) is replaced by 

M 

f/= exp(-i0™i/fc„), (5) 

m— 1 

where dm — At/hum and km labels the Hamiltonian 
from set S that is to be used for the mth pulse. As the 
system is stationary in the absence of laser interactions, 
the system Hamiltonian i/sys was omitted in Eq. 

Using the gradient-ascent method for updating the 
control amplitudes Om increases the performance func- 
tion but leaves the pulse ordering defined via the indices 
km unchanged. Thus, for the optimum pulse order to be 
included in the configuration space, the number of pulses 
M needs to be much larger than the expected minimum 
number of pulses finally realizing the target operation 
CAargct- Therefore, the search algorithm has to be com- 
plemented by a penalty function like <I>p = X]m=i l^'mPi 
< 7 < 1, that tries to eliminate short pulses that 



do not contribute much to increasing the performance 
function $. The functional form of <I>p assures that a 
change in the length 9m, of the m'th pulse by d6m pe- 
nalizes already short pulses much more than longer ones 
since d^p = sign{9m)l\0mP^^ dOm- In the optimization 
routine used for finding the pulse decompositions pre- 
sented in section IIVI the exponent 7 ranged from 0.5 
to 0.8. Now, the performance function $ is replaced by 
$ = $ — a$p where a is a suitably chosen weight. 

For updating the pulse lengths 6m, we do not calculate 
the gradient of $ and move in the direction of steepest 
ascent but perform consecutive one-dimensional maxi- 
mizations of 9m instead. This has the advantage that 
the step size e of Eq. ^ can be individually adjusted for 
the different directions by considering also the curvatures 
d^'t/d9fj^. For a negative curvature, the pulse length 9m 
is updated by going to the maximum of the parabola 
approximating the $ whereas for a positive curvature a 
fixed step size is used for updating 9m- 

To avoid becoming trapped in a local maximum of the 
performance function, we combine the uphill search algo- 
rithm with elements of simulated annealing. Instead of 
choosing the pulse length 9m corresponding to the maxi- 
mum of the parabolic approximation to $, the algorithm 
samples the region around the maximum by randomly 
choosing a pulse length 9m = 9^ + A9 where A9 is ran- 
domly drawn from a normal distribution with probability 
density cx exp{{d^^/d9l^){A9)'^ /Tes) where the effective 
temperature T^s determines the spread of the distribu- 
tion around 9^. In the course of the optimization. Toff 
is lowered to zero. In addition, the algorithm tries to in- 
troduce new pulses into the sequence from time to time 
to achieve a variation of the pulse order. Unless other- 
wise mentioned the program is started from a random 
sequence of pulses of sufficient length. 

The computational overhead for performing an up- 
date of the pulse amplitudes is the same for a method 
following the gradient and for the M consecutive one- 
dimensional optimizations of the pulse lengths. Choosing 
the one-dimensional optimizations allows us to calculate 
the curvatures 8^^/89^ at no additional cost as well as 
to include the annealing technique. It would be prefer- 
able to use all the information encoded in the elements 
d^'t /d9md9n of the Hesse matrix, however, we judged the 
computational cost amounting to 0{M^) matrix multi- 
plications to be prohibitively high. 



IV. EXAMPLES 

The optimization routine was used to search for decom- 
positions of unitary transformations of interest in systems 
of three to five qubits. In the following subsections, ex- 
amples will be given for pulse sequences found by the 
program. Interestingly, it turns out that in most cases 
the sequences consist of pulses having pulse lengths that 
are simple fractions of tt even though the optimization 
routine was allowed to vary the pulse lengths continu- 
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ously. Moreover, the pulse sequences listed below realize 
the desired target operation not only approximately but 
exactly. 

Depending on the initial pulse sequence and the values 
of the parameters controlling the optimization process, 
the optimization algorithm can converge to different so- 
lutions realizing a target operation J/targct- This demon- 
strates that we cannot be sure that the pulse sequences 
found by the program necessarily represent the optimum 
solution. However, as we are interested in discovering se- 
quences of practical interest to be used in experiments, 
this is hardly a drawback. 

In the following, we simplify our notation by using 
the short-hand notation ui''\e) ^ [0]^, Ux{0) ^ [e]x, 
^Msi^) ^ [(^]xx to achieve a convenient and compact 
representation of the pulse sequences. Pulses are sepa- 
rated by hyphens and temporally ordered from left to 
right. Whenever a sequence is given, the number of 
qubits it is operating on is mentioned in the text. 



B. CNOT gates on more than three qubits 

In general, a sequence realizing a two-qubit gate in a 
system with = 3 qubits will not correctly function in 
the case iV > 3 as it risks to entangle the spectator qubits 
with each other. On the other hand, a sequence realizing 
a two-qubit gate for N = 4 which does not contain any 
phase shift gates on the spectator qubits is also applicable 
to > 4 because of the symmetry of the interactions 
between the different spectator qubits. A pulse sequence 
realizing a CNOT gate for A^ > 4 is given by 



[fix 



l^lxx - 



[^\X - [TTJ^ - [--\XX 



(6) 



where again the qubit 1 controls the target qubit 2. 



A. CNOT gates on three qubits 



C. Further multi-qubit operations 



In a system of three qubits, a CNOT gate between two 
qubits is realized by the sequence 



i^lx 



l^lxx 



l^\x 



l^\x 



\xx 



where the target qubit 2 is controlled by qubit 1, cor- 
responding to the operation U — {I + a 



(1) 



J2) 



A unitary transformation consisting of two CNOT op- 
erations with qubit 1 controlling the other two qubits 
is described by the unitary operation U = {X ^ a^z^ + 

(T^x^ c^x'^ — o'i^''(Ti^-'cri^-')/2. One way of decomposing it 
into elementary operations is given by the sequence 



-j\xx 



\xx 



V^\x. 



\x 



\x 



by 
of 



As an example of a pulse sequence found 
the search algorithm that is not composed 
pulses having pulse angles which are simple ra- 
tional fractions of tt, we present another decom- 
position of the threc-qubit operation [/ given by 



■2 XX 



'2 XX 



"1 



[/3 



1 XX 



«1 



For ai w 0.712l7r, w -0.424l7r, /3i « -0.212l7r and 
/32 ~ 0.35607r, it also realizes up to an unimportant global 
phase two CNOT gates on three qubits with qubit 1 con- 
trolling the other two. The program having provided the 
pulse angles, we found in a second step that the angles 
satisfy the algebraic relations a\ = 2/^2, a2 = ir — 4/^2, 
/3i = 7r/2 - 2/32 with = In - \ arcsin(\/5 - 2). 



Apart from CNOT gate operations, we also searched 
for a decomposition of a quantum Toffoli gate operation. 
The sequence HI] 



'2' 



Jxx 



IX 



2' 



JXX 



Jxx 



Jx 



is applicable to a system of three qubits and flips the 
state of the third qubit depending on the state of qubits 



(3X 



.(3) 



1 and 2, thus realizing the operation U 
(aW+af -^i;VF)(I-a(^'))/4. 

In a system with an even number N = 2M of qubits, 
a mapping between the Bell basis and the product state 
basis for each pair of qubits (2m— 1, 2m), m = 1, . . . M, is 
of interest as it could be used for measuring the multipar- 
tite concurrence of an M-qubit quantum state available 
in two copies [23l | and for implementing entanglement 
purification protocols. For the case of four qubits, the 
sequence 



[4JXX 



Ml Ml - [^]xx 



Ml - Ml 



realizes the desired mapping described by C/(i_2 3-4) = 

exp(-zf(ai^Vi^)+ai^Vi^))). 

The approach can be extended to higher numbers of 
qubits. Executing the sequence 



Jxx 



Ml - Ml 



[ttJxx - m. 



twice, realizes the unitary transformation C/(i_2.3_4, 



exp(-if Em=i 



(2n 



-1) (2m)s 



(1 

for the case A^ 



6,8. 
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spin flip errors 



spin flip errors 



a|0>+|3|1>- 



io>-e 

|0> 



encoding 



-y h- 



|0> 
|0> 



-V A- 



lo>- 



. ancillas 



|0>- 



syndrome 
detection 



correction 
step 



= U 



QEC 



FIG. 1; Repetitive error correction of spin flips using three 
qubit for encoding the logical qubit and two ancilla qubits for 
error syndrome detection and coherent correction of errors. 
The optimal control algorithm looked for a decomposition of 
the unitary operation enclosed in the dashed box comprising 
syndrome detection and error correction. 



Moreover, this approach could also be used to create lin- 
ear cluster states and to realize the Hamiltonian of a 1-D 
Ising model by concatenating the unitaries C/(i-2.3-4....) 
and C/(2_3,4-5,...)- 



D. Quantum error correction 

In experimental quantum information processing, first 
steps have been taken to demonstrate quantum error cor- 
rection [13, HI] using three-qubit codes. A major step 
forward would be the realization of repetitive error cor- 
rection where the quantum information remains encoded 
in a logical qubit all the time. In the example shown in 
Fig. [ll a qubit state a|0) -|- /3|1) is encoded in a logical 
qubit consisting of three qubits as I^'l) = a|000)-|-/3|lll). 



The circuit detects single bit flip errors Ox 



(m) 



1,2,3, 



by means of two additional ancilla qubits for syndrome 
detection and corrects the errors coherently by applica- 
tion of quantum Toffoli gates in order to restore the state 



Using gate decompositions similar to Eq. ([6]) for con- 
structing the unitary transformation J7qec that detects 
and corrects errors from CNOT and Toffoli gates would 
result in a pulse sequence with more than 100 pulses that 
does not seem to be realizable with current technology. 
Therefore, the only practical approach seems to be to 
search for a gate decomposition of the complete opera- 
tion C/qec- However, as the time needed for finding gate 
decompositions scales exponentially with the number of 
qubits, the task of decomposing the five-qubit operation 
[/qec might seem quite challenging at a first glance. But 
actually, there is a whole class of operations J7qec equiv- 
alent to ?7qec accomplishing the error correction pro- 
tocol. To see this, note that we only require C/qec to 
perform the mapping [/("'^ |00) |V'l)|V'a"^) where 
C/(™) e {J,cri^\(Ti^\(Ti^^} and are arbitrary or- 

thonormal vectors describing the ancilla state at the end 
of the correction step. Therefore, any valid transforma- 
tion C/qec can be put into the form 



C/qec -foo 



E„(|000)(000|C/™®|V'Xo)(00| 
+ |lll>(lll|C/„® |7^Xi)(ooi) 



(7) 



where the ancilla states needs to satisfy the constraint 
= ^/,™^. Here, Pp^ = 1® |00)(00| is a projector onto 
the initial ancilla state. The condition (?/'a,o|'0a,i) — 1 
expressing the constraint ^/'™q = assures that no 

information about the state of the logical qubit can be 
obtained from detecting the ancilla state. To search for 
a gate decomposition, the performance function can now 
be modified by replacing the trace Tr(C/^C/QEc) by <& = 

Re(E„.(V'!roVS». 

Using this approach, the search program found the fol- 
lowing pulse decomposition C/qec for realizing an opera- 
tion equivalent to C/qec : 



'2' 

rTT, 



^1 



\x 



^1 



[ttIxx 



\xx 



\xx- 



^-.\xx 



Stt, 



Stt, 



\YY - \-j\X - V7^\ 

' \-l\xx - \Al - 



Here, the uphill search was extremely slow when starting 
the search algorithm from a random sequence of pulse. 
However, starting the algorithm from a sequence con- 
structed from gate decompositions of the CNOT and 
Toffoli gates and initially driving it away from this unde- 
sired solution by increasing the effective temperature Tog 
proved to be an effective strategy for finding an improved 



solution. 

V. SUMMARY AND OUTLOOK 

In conclusion, we have developed an algorithm based 
on optimal control techniques for finding decompositions 
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of unitary transformations into finite sequences of pulses 
that correspond to the apphcation of Hamiltonians drawn 
from a given set. For a set of Hamiltonians that are of 
interest for ion trap quantum computing, the algorithm 
provides gate decompositions that would be difficult to 
find otherwise. 

Looking for gate decompositions involving parallel two- 
qubit interactions on a small number of ions provided in 
many cases (as for example for the quantum Toffoli gate) 
a pulse sequence much shorter than what would have 
been possible by using sequential two-qubit interactions 
and single qubit gates. In a few other cases, however, 
this strategy did not pay off (as we noted when looking 
for a decomposition of the quantum Fredkin gate). 

While our investigation was limited to a particular 
basis set of Hamiltonians S, the program could be 
adapted to search for gate decompositions using other 
sets that might be more relevant for realizations of 
quantum computing in other physical systems where, for 
example, the natural entangling interaction is given by a 
ViSWAP gate ^ or an exchange interaction W^. For 
example, in the context of ion trap quantum computing. 



another interesting question is whether the algorithm 
could be modified to make it applicable to the case 
of quantum gates realized by a tightly focussed laser 
interacting with a single ion at a time. In this case, the 
Hilbert space would be comprised not only of the qubit 
states but also of the harmonic oscillator the qubits 
are coupling to, a configuration not only found in ion 
trap quantum computing but also in cavity-QED setups 
where an electromagnetic field mode is interacting with 
cold atoms or superconducting flux qubits. 
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